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I. INTRODUCTION 


There are at least three components which contribute to the total energy density of 
the Universe. These are non-relativistic matter, relativistic matter and dark energy. Dark 
energy, which constitutes about 70% of the present energy density of the Universe is one of 
the most discussed and hot topics of this era. Several possible candidates (explanations) of 
inflation and dark energy are discussed in the literature, a few of them are: the cosmological 
constant (A), for a brief review see PQ |2|; modified gravity, for a good discussion see 011!; 
scalar held models (e.g. quintessence, K-essence, tachyon held, phantom (ghost) held [SHT], 
dilatonic dark energy, Chaplygin gas) [2] and vector helds|H| |9]. 

Scalar, vector and tensorial helds came into existence with the development of quan¬ 
tum held theory and were hrst used by Alan Guth [10], Andre Linde [11] (plus references 
therein for discussion), to propose inhation as a possible solution to the horizon and hatness 
problems. For a review on inhation see [12j . 

The homogeneous, isotropic Universe model is described by the Friedmann-Robertson- 
Walker (FRW) metric and its dynamics is described by the Friedmann equations which are 
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where a is the scale factor, if = a/a is the Hubble parameter, p is the total energy density, 
K is spatial curvature and P is the pressure. Where, cUefr = ^de ^de + 
and the barotropic equation of state P = ujp has been used. In the equations, natural units 
h = c = 1 and (SvrG)”^/^ = Mp = 1 have been in above and subsequent equations. From 
eq. (|^, we see that a > 0 when cuefr < 

Using the law of conservation of energy, it can easily be derived that 
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From eq. ([^ scalar helds are categorized into two categories on the basis of their equation 
of state (EoS) parameter ca: 
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The Quintessence field: These fields lie in the domain —1 < a; < —1/3. The main 
feature of this field is that its energy density decrease as in relativistic (ca = 1/3) or 
non-relativistic matter {oj = 0), but the field could be a possible candidate for dark 
energy since oj < —1/3 [2]. 

The Phantom field: These fields he in the domain a; < —1 and have the property that 
their energy density increases as the Universe expands and the scale factor ‘a’ becomes 
infinite in a finite time. These type of fields have negative kinetic energy, which allows 
u to take value less than -1 [2]. 

When a; = — 1 the energy density does not change as the Universe expands and remains 
constant. The deviation of u from —1 tells how us energy density, p, of the field changes as 
the Universe expands. 

The scalar fields used to explain accelerated expansion of the Universe are taken to be 
new fields (with no connection with Particle Physics). However, the introduction of new 
fields with no experimental basis other than “explaining” one observed phenomenon is too 
reminiscent of epicycles or the per-relativistic aether. Only if one can exclude fields contained 
in the standard model, or a minor extension of it, would it be justifiable to introduce such 
exotic proposals as phantom or quintessence fields, which I argue is necessary and essential 
to made a unified model of matter and energy. 

In this paper, I assume that the present state of the Universe is described by an inert 
Higgs doublet model. If the phase transition in the second SU(2) doublet, (^ 2 , (which has 
zero vacuum expectation value (vev)) has not occurred yet, then the component fields of the 
doublet can be considered as dark energy candidates. In the case when the present Universe 
is described by an inert Higgs doublet model after the second phase transitions different 
forces and particles would arise (depending upon the VeV of the second doublet) which can 
be avoided by taking the parameter 77122 Higgs potential zero thus second phase transition 
never occurs. The results obtained are presented which favor that accelerated expansion 
of the Universe is quite possible (in quintessence regime) [T^. Previously, E. Greenwood et 
al in na considered the appearance of a second vacuum in the SM Higgs potential. They 
showed that if the second order electroweak phase transition is followed by a first order phase 
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transition (which they were able to achieve for a certain range of parameters of their model) 
then this phase transition can canse the late time acceleration of the Universe. Carroll et al 
have also given an independent phantom field model in [5] where they have also discnssed the 
stability of energy density along with phantom field as a possible dark energy candidate.V. 
K. Onemli and R. P. Woodard have also showed in that Qnantnm effects in Cosmology 
conld canse the violation of the nnll and weak energy conditions (i.e. p + P < 0 in their 
model) cansing a; < — 1 on cosmological scales withont introdncing any ghost, phantom, etc.. 

The paper is organized as follows; in section 2 we review the inert donblet model and 
obtain constraints on the inert donblet model’s potential to be positive along with the calcn- 
lations of the minima of the potential. In the next section phase transitions of inert donblet 
model has been discnssed. Section 4 is devoted to obtaining the answer to the qnestion; 
whether Higgs field can be a possible candidate for dark energy in inert donblet model? In 
section 5, we discnss the decays of Higgs field(s). Section 6 is where we conclnde the paper 
with fntnre work possibilities. 


II. THE INERT DOUBLET MODEL 

The electroweak symmetry in the standard model (SM) of Particle Physics is broken 
spontaneonsly by the non-zero vev of the Higgs field(s) via Higgs mechanism. The Lagrangian 
which describes any model in Particle Physics is 

^ + ^mggs • (4) 

Here, .^g/^ is the SUc{S)<^SUl{2)<^Uy{1) is the SM interaction of the fermions and gange 
bosons (force carriers) [15] , ^y is the Yukawa interaction of fermions with the Higgs 
field(s)[T6] and ^Higgs is the Higgs field Lagrangian where 

■^Higgs = Th — Vh , (5) 
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Th being the kinetic term of the Higgs £eld(s) and Vh being potential of the Higgs £eld(s) [IT] . 
Here, 

Th = 

( 6 ) 

+ [x(-Dl/i0l)^(-D2^02) + X*(-D2^02 )^(-Di^</>i)] , 

and 


Vh= V 1 + V 2 + V, 


int 


(7) 


= Pi exp(Aimii0j0i) + p 2 exp(A2m220202) + Ps exp ^-A3Ai(0l0i)^ 

+ P4exp +mi2(0i02) +m?2*(020l) + -^3(0101) (0202) 

+ '^ 4 ( 0102 ) ( 020 l) + X ['^5(0102)^ + -^ 5 ( 020 i)^] + -^ 6 ( 010 l) ( 0 l 02 ) 


( 8 ) 


+ A6(010i)(020i) + A7(0202)(0l02) + A;(0^02) (0^0l) , 
where Vx and V 2 in eq. ([^ are the Lagrangian of Higgs held 0i (given by hrst two terms of 
RHS of eq. (|^) and 02 (given by 3’’'^ and 4*^ terms of RHS of eq. (|^) respectively, Vint is 
the interaction Lagrangian of fields 0i and 02, given by remaining terms. 


D2tL — dfiT + ^^Bfj. , 


A _ 


Vi - i^Xi + 


Vi + >'Xi + ^i 
The dimensions of different parameters are 

[pr^ = [A.] = [A]^ [ml] = [L]-\ [0,] = [L]-^ and [A,] = [L]0 

Here “L” denotes the length. The helds 0^, 0“, p* and Xi the hermitian Higgs helds (0^ 
are charged whereas other helds are neutral), Ui is the vev of the doublet 0*. 

The Higgs vacuum energy of the potential given by eq. ([8| is 


(9) 


Evac = Pi exp +p2exp Qm|2A2^'|^ Tpsexp QaiAbz/^ 

+P4exp ^-A2A4i/|^ + -Asi/fz/l + -A4i2^z/| + - (A5 + Ag) 

+ ^ (-^6 + -^e) ^ 1^2 + ^ (-^7 + -^ 7 ) ^1^2- 
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We would be choosing z /2 
doublet model, then 


0 here (explained in section IIB), which corresponds to the inert 


p2 + Pa + Pi exp 




+ P 3 exp 
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Since we are looking to provide the accelerated expansion to the Universe from Higgs 
£eld(s) we want the Higgs field(s) to live longer than the age of the Universe (or at least not 
less than the present age of the Universe). Thus, stable Higgs £eld(s) is(are) required for the 
scalar held dark energy model. This can be achieved by imposing a discrete Z 2 symmetry 
0 —)■ —0. When Z 2 symmetry is broken by (0j0j) type terms then it is said to be softly 
broken and when Z 2 symmetry is broken by type terms then symmetry is said 

to be hardly broken. The terms containing m ^2 describe the soft symmetry breaking of Z 2 
symmetry, the terms Xq and A 7 describe the hard symmetry breaking of Z 2 symmetry. In 
the absence of these terms along with no cross kinetic term i.e. x = 0, the 2HDM’s Higgs 
Lagrangian has a perfect Z 2 symmetry [18] . There are two Z 2 symmetries corresponding to 
the doublets: 


I: (j)i —)■ - 01 , 02 —t 02 , (11) 

H: 01 -)■ 01 , 02 -)• -02 • (12) 

For Z 2 symmetry see [T 8 ] . 

It should also be noted that even if the Higgs £eld(s) has(have) decayed off its vacuum 
energy given by eq. (|^ would still be able to provide accelerated expansion provided it 
dominates at some time in the life span of the Universe. 


A. Positivity constraint on Higgs potential 

The stability of the vacuum given by eq. ([^ implies that the potential must be positive 
for all asymptotically large values of the Higgs held. Even though Ui and V 2 remains positive 
for all values of helds considering pi > 0, we should End out the constraints on the parametric 
product Aim^i, A 2 m 22 , ^AsAi and IA 4 A 2 when the exponentials are approximated upto a 
hnite order (4th here). The conditions obtained should match on to the conditions for the 
general 2HDM potential given in [T 8 ] . 
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Some conditions can be obtained to achieve the stability of the vacuum by introducing 
new parameters as mm, 


|0i|=rcos7, 1021 =?" sin 7 , 


401 


= pexp{l9) , 


I0l| I02| 

where 7 G [ 0 , 7 r/ 2 ], p G [0,27r). Using these transformations in the potential within its 
expanded form upto forth order (with > 0 and hence will be omitted in eq. (HI) , the 
potential becomes 

U = cos^ 7 + ^2 sin"^ 7 + A 3 cos^ 7 sin^ 7 + X^p^ cos^ 7 sin^ 7 

+ Asp^ cos^ 7 sin^ 7 cos 26 + 2 pcos 7 sin 7 cos 6 '[A 6 cos ^7 + A 7 sin^ 7 ] 


(13) 


the positivity of eq. (13) implies that 

Ai > 0 , A2 > 0 , A3 > —2 A^A2 , 
when Ae = A 7 = 0 , then we also have 

A 3 + min[ 0 , A 4 - lAsI] > -2 ^A'^Ag , 

where A) = ^ (P 3 A 3 A 1 + pi(Aim? 02 ) and ^2 = ^ (P 4 A 4 A 2 + p 2 {^ 2 ml 2 f). 


(14) 


(15) 


B. Minimizing the Higgs potential 


The extrema of the potential are found by taking 
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</'2 = (< A 2 > 
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and 
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</'i = (0l> 
</'2 = ( 02 > 


The most general solution of the conditions (16) is 


dVH 

502 


</'2 = (</>2> 


0 . 


(16) 


(0i) = 


1 

71 



and 



The hrst solution of extrema has been taken to be similar to the Higgs vacuum in SM and 
the second one is the most general which can occur. One needs to keep in mind that now 
|z/|| +v?. In the electroweak symmetry breaking situation where (0i) 7 ^ 0, choosing 
the z-axis in weak isospin space such that i^i > 0 . 
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When u ^ 0, the non-zero value of u will contribute to the “charged” type dark energy, 
which has not been observed. To avoid this, we would take u = 0. From the extrema 


conditions given by eq. (16), we can determine the values of ui and z /2 psusp, solving eq. 
(16) for the potential given by eq. ([^ leads to 

/Aim'll 




exp j exp j - 2 Re(m? 2 )^ 


T (-^3 + A 4 -|- Re(A 5 ))z^| -|- 3Re(A6)z/4Z/2 -|- Re(A 7 ) — 

z/i. 

p4A4A2r'|exp -2Re(m?2)^ 

T (-^3 + A 4 -|- Re(A 5 ))z^^ -|- Re(A 6 )— + 3Re(A7)z/4Z/2 

Z ^2 


3n 


Z ^2 


2 p 2 A 2 ml 2 exp 


V 2 




= 0 , 


= 0 . 


The solution of the above equations for ui and i >2 is practically impossible since the expo¬ 
nentials contain Ui and z/ 2 . The vev’s are approximated by truncating the potential eq. (|^ 
upto forth order. This is also a good approximation in the sense that the results of u's would 
also be true for the general 2 HDM potential vev’s of [IH]. Applying minimization conditions 


given by eq. (16) we get. 


z^i 


2 piAim^;^-|-2 A'^z/^ — 2 Re(m^2)~ + (A3 + A4-|-Re(A 5 ))z/, 


Z ^2 


z^i 


2p2-A-2'zri22 T 2A2Z/I — 2Re(?Tr^2)— T (A3 + A4 -|- Re(A5))z/^ 

Z ^2 


3Re(A6)z/iz/2 -I- Re(A7) — 

z^i 


(17) 


= 0 , 


z/r 


-|- Re(A 6 )-h 3Re(A7)z/iz/2 

Z ^2 


(18) 


= 0 , 


where A) = ^ (P 3 A 3 A 1 -h pi^AimhY) and A 2 = ^ (P 4 A 4 A 2 + p 2 (A 2 m| 2 )^) as given before. 

Since Z 2 symmetry is very important in determining the phenomenology of the theory 
and with exact Z 2 symmetry the lightest Higgs boson (held) is stable. In this article, the 
potential which is Z 2 symmetric will also be taken and the possibility of any Higgs held(s) 
to be the candidate for the dark energy will be examined, if possible it would give an ever 
expanding Universe while accelerating. Thus, imposing the Z 2 symmetry [18] we have. 


X = 0 , 


77142 = 0 , 


Ae — 0 , 


A7 — 0 , 


(19) 

















Using eqs. (17, 18, 19), we get four solutions for vi and 1 ^ 2 , which are 


ul = 0 , 


K = 0 , 


= 


u, = 


a; 

2{2piAiX'2Tnl^ - p2A2\‘ii^ml^) 


4 A]^A 2 A345 

where A 345 = A 3 + A 4 + Re(A 5 ). 


z/| = 0 , 


Z/n = 


p2A2ml2 


a; 


z/| = 0 , 


Z/n = 


2{2p2A2X'iml2 - piAiXsi^ml^) 


4 A 4 A 2 A| 


(20) 

( 21 ) 

( 22 ) 

(23) 


345 


Electroweak symmetric vacuum This vacuum solution is given by eq. (20) represents 


the electroweak symmetric extremum, in this case the complete Lagrangian respects Z 2 
symmetry with respect to both fields given by eq. ( [II| and eq. (12) i.e. I and II, here 
mh > 0 and m 22 > 0. The doublets are 



Vi + IXi 


and 



Inert vacuum-I Eq. (21) will be called inert vacuum-I. The doublets in this vacuum are 


Vi + 


and 


V2 + 1X2 + ^2 


in this type of vacuum solution, the fields and X 2 are the Goldstone bosons and all other 
fields are physical. In this case, the complete Lagrangian respects Z 2 symmetry given by eq. 
( pT| and only the Higgs Lagrangian (not the complete Lagrangian because of the Yukawa 
Lagrangian since 02 is not contained in Yukawa Lagrangian) violates given by eq. (12). 


Inert vacuum-II The doublets in this type of vacuum, given by eq. (22), are 


01 = 


5 4>2 — 



Vi + i^Xi + G 


V2 + 1X2 


(24) 


the fields 0^ and xi ^-re the Goldstone bosons whereas other helds are physical. The Yukawa 
interaction is described by the interaction of held 0 i with fermions (02 does not couple to 
fermions but it can appear in loop process). 
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The Higgs and Yukawa Lagrangian violates Z 2 symmetry given by eq. 0 whereas 


respects given by eq. (12). Thus, parity in held 02 is conserved this makes the lightest 


particle/held of doublet 2 stable. If there is no phase transitions (which brings Universe 
in the mixed vacuum state) in the future of the Universe or before the phase transitions 
in doublet 02 then this lightest particle is a good candidate for dark energy, giving rise to 
an ever accelerating expanding Universe. Without the stability, we would have a model for 
accelerated expanding Universe for a limited time depending upon the decay width (life time) 
of the dark energy £eld(s), because if £eld(s) can decay then it ceases to provide accelerated 
expansion after its decay to the other helds only if its vacuum energy has not dominates 
before its decay. 

The masses of the helds in this vacuum are 


< = 2 A>? , 


ml = p 2 A 2 m |2 + 


= P2^2mi2 + 

= P2^2ml2 + ^ 


A3 + A4 + Re(A5) 

2 ^ 

A3 + A4 — Re(A5) 


'1 ) 


1 5 


(25) 


where vi = 


</2Gi^ 


246GeV. 


In this paper, I am interested in this type of vacuum. 

Mixed vacuum After symmetry is broken, the Lagrangian violates the Z 2 symmetry 
given by eqs. ( [TT| and (12). The physical helds in this model are the combination of helds 
from 01 and 02. The complete theory and phenomenology of this model is discussed in [22] . 


Mixed vacuum solution is given by eq. (23). 


III. PHASE TRANSITIONS IN INERT DOUBLET MODEL/TWO HIGGS DOU¬ 
BLET MODEL 

When we look for dark energy to be some physical held(s), then it becomes essential to 
also look for its evolution in the history (cooling) of the Universe. 

According to [21], in the quantum held theory at non-zero temperature the terms of 


10 









the quadratic terms evolve with temperature as 


where 


ml, + -c,T\ 

2 . ^ rr2 

^22 + 2 *^^^ ’ 


_ 6Ai + 2A3 + A4 ^ ?)gl + g'^ ^ gl + gl 

S + ~^3r~ + “ 8 “ 


_ 6 A 2 + 2 A 3 + A 4 ^ + g2 

- 12 + 

here A'^ = ^ (P 3 A 3 A 1 + pi(AimiA 2 = ^ (P 4 A 4 A 2 + P 2 (A 2 ?Ti 22 )^) as above and pi ,2 and 
g[ 2 are the Electroweak gauge couplings with doublet and 02 respectively, g, ^ 1 and 
gb ~ 0.03 are the top and bottom quark Yukawa couplings with the doublet 0i in inert 
doublet model respectively. 

In general, ci and C 2 can have any sign but the potential positivity implies that (in any 
situation) 

Cl + C 2 > 0 , 

with the above mentioned positivity constraints, the case when ci > 0 and C 2 > 0 will be 
considered. 

The Higgs potential with the new quadratic terms now become 

Vh(.<Pu<P2,t) = K„o + yi(r"-r=;)(yi.#.,) + 02(r2-r,^)(4.^2) 

+ A'i(0l0l)^ + A2(0202)^ + m?2(0102) + m?2*(020l) 

+ A3(0i0i)(0202) + A4(0|02)(020i) + 2['^5(0102)^ + 

+ (010l) [-^6(0102) + K>i.4>2(l>l)] + (0202)[A7(0l02) + A7(020l)] 

+higher order terms , 


In this case corresponding to the two different vev’s of the helds 0i and 02, there would be 


2 piAimfi 


occurring two different critical temperatures Tc, = W-——-—^ and Tc^ = \l - ——^—— 

V Cl V C2 

respectively. In this analysis we require that the phase transitions in doublet 02 occur after 
phase transitions in 0i, thus Tc^ > Tc^- 


2 p2A2m22 
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There are three possible situations depending upon the value of temperature, 


T >T, 


Cl" 


For T > Tc^ the symmetry is not broken in either field and the phase transitions has 
not occurred. Thus, in this situation minima of both helds lies at zero given by eq. 


(20). All fields are massless in this temperature range i.e.T > T^. Since, symmetry is 
not broken in this limit all the fermions and bosons are also massless. 


<r<T,,: 

This is the most interesting situation as the symmetry is retained in field 02 and is 
broken in 0i. In this situation due to phase transitions in 0i all helds (fermions and 
all bosons) gets masses, but since the phase transitions in held 02 has not occurred so 
it still is a possible candidate for dark energy held. This situation is also interesting 
because this article assumes that the present observed Universe happens to occur in 
this domain. The expansion of 0i around its minimum vi leads to exactly the same 
phenomenology as that of the SM. 

T<T,,: 

When T < symmetry has got broken spontaneously in both helds 0i and 02 at 
temperatures T^i and Tc 2 respectively. In this situation, both helds gets perturbed 
around their true minima. This perturbation give masses to all the mixed helds of the 
doublets 01 and 02 (except the CP odd and charged held of angle 0 = tan“^ ^ rotated 
new Higgs doublet Hi which are the Goldstone bosons). This should be expected to 
occur some time in the future considering m \2 ^ 0. 


Constraint on parameters from phase transitions 

If we require that the phase transitions in the doublet 0i (which is acting now to be SM 
Higgs doublet) occur at the same temperature as in SM Higgs doublet, we can constrain the 
parameters of the potential. Note that this condition is not necessary but an interesting one 
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which can probe this model (one can also relax this condition). The critical temperature at 
which phase transitions occur in SM is [ 23 ], 


rTi2 _ 

C 


4 A 1/2 




where A = 0.1305 is the quartic coupling of Higgs doublet in SM, g = 0.6376 is the Higgs 
doublet coupling with SU(2) gauge group and g' = 0.3441 is the Higgs doublet coupling with 
U(l) gauge group. With the values of parameters given above we get = 230.3186GeV. 
Setting Tci = gives, 


-2 Pi Ai mh 


4 A 1/2 


6 Ai + 2 A 3 + A 4 3gj + g'l gj + gl \ 2 A+-o 2 + io'2 
\ 12 ^32 ^8 ^44 


taking square and using equation ( 22 ) in, 

2 a) z/2 


4 A z/2 


6A1 + 2A3 + A4 3 g >2 + g2 g2 ^ ^ 

12 32 8 4 4 

here i/2 = i/2 = 246 GeV2, Aj’s are the parameters of inert doublet model potential, gi 
and g'l are the couplings of doublet 0i with the SU(2) and U(l) gauge group respectively. 
Simplification of the above equation with the numerical values lead to 


2 A 3 T A 4 — 0.7869 , 


Here, in deriving the eq. 
have been used. 


(29), Ai 


0.1305, gi 


0.6376, g'l 


0.3441, gt 


1 and gb 


(29) 

0.03 


IV. HIGGS FIELD(S) AS DARK ENERGY FIELD(S) 

The second Friedmann eq. (|^ tells that accelerated expansion will occur when Weff < — | • 
For the field (^2 to be the dark energy field, it must bring cUes < — | in the history of Universe 
(in fact just now Z 0.37 with Gm = 0.3 and Gde = 0.7) which tells that the symmetry 
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in this field should not be broken. For this purpose, we need to solve the Euler Lagrange 
equations in some background, which are 


Q f ^iV 9-^Higgs) \ _ d{\/ 9-^Higgs) _ ^ 

^ V ) dij. 


where t/ji are different fields of doublets 0i and 02 - 

The solution of Euler Lagrange equations of motion in FRW Universe {y/—g = ct{t)^) for the 
fields 02 ) V 2 and X 2 given above is 


V2 


3-^72 + - 
a 2 


V 2 { (As + A 4 + A 5 ) z/2 + 2 ml^A2p2 exp ( -ml 2^2 (xl + hi + 20f) 


1 


+A 2 A 4 P 4 xi exp ( AA 2 A 4 (X 2 + hi + 20 f ) + 2 A 2 A 4 p 40 f exp ( AA 2 A 4 [xl + hi + 20 f) 


. 2^2 


2\2 


+A2A4p4hi exp -A 2 A 4 (xi + hi + 20f) 


= 0 


(31) 


:A 2 + 3 ^x 2 + ^ (x 2 (^(As + A 4 - As) z /2 + 2 ml 2 K 2 p 2 exp Qmi 2 A 2 (xi + hi + 20 f) 
+A2A4P4?7| exp ( ^A2A4 (xi + hi + 20 f + 2X2A4p4<pf exp ( ^A2A4 (xi + hi + 20 f 


+A2A4p4xi exp ( -A 2 A 4 (xi + hi + 20f )■ 


= 0 


(32) 


0 ^ + 3 ^ 0 ^ + ^ I^Asz/^ + 2 ml 2 A 2 p 2 exp Qmi2A2 (xi + hi + 20f ) 
+A2A4P4 (xi + hi + 20 f) exp f JA 2 A 4 (xi + hi + 20 f 02 = 0 


(33) 


To avoid the confusion between positively and negatively charged Higgs boson, ‘c’ instead 
of + or — has been used. 

The energy density and pressure after expansion of 2HDM Higgs Lagrangian for physical 


14 



fields is 


pHiggs/Pniggs — 


+ if ± 


7 (As + A4 + As) V ^2 P 1; (A3 + A4 — As) z/ X2 


+ ^A3Z/Vf + P 2 exp (^ml^A 2 xl + ^m| 2 A 2??2 + ml2^24>f 
+P4exp QA 2 A 4 X 2 + ^A2A4r7| + ^A2A40^‘^ + ^\2A^xh2 
+ -j^^2^iX2^2^ + 2^2^4V2^2^ 


( 34 ) 


For the cosmological evolution of the helds 172, X2 and 02, the above equation of motion are 
solved with the Friedmann equations numerically in the flat Universe (k = 0). The initial 
conditions used are 772 ^™ = ATp, X2ini = Mp, 0|.„. = 0, i) 2 ini = 0, X 2 ini = 0 and = 0. 

Note that after imposing the Z 2 symmetry there are hve parameters which determine the 
masses of Higgs helds. Withm^^ = musM — 125.7GeV and constraint from phase transitions, 
we have only two equations to determine parameters values[21]. The Higgs helds masses in 


this analysis were calculated using eq. (25). After imposing = musM — 125.7GeV 


and constraint from phase transitions eq. (29) some the arbitrary choice of parameters was 
taken [25] since there were only two eqnations to determine hve nnknowns. 

Masses of Higgs bosons in the analysis are taken to be 


= 6.925 X 10“®^GeV, = 6.925 x 10“®^GeV, m^± = 154.305GeV, 

jnst to set the evolntion of the energy densities as observed. Bnt too small masses of lightest 
Higgs helds were reqnired to indnce high damping in the held’s oscillations which is reqnired 
to ensnre that cjefr get less than —1/3 only once on the cosmological scale and whenever it 
gets it never gets valne greater than —1/3. Since the vacnum energy of the Higgs is too 
small (!l(10“^^^i7p/Lp), it asks for very small values of parameters to achieve over-damped 
oscillations. Thus, the hne tuning problem related to the smallness of held mass in scalar 
held models also exists in our model too. 

To make the Universe evolving in the vacnnm states given by Inert vacunm H, m \2 = 0 
has been chosen which implies Tcj = 0. In this sitnation, phase transitions in the second 
doublet will never occur (no appearance of odd forces and particles would be observed) and 
the Universe will keep on accelerating forever. 
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The solution of the eqs. (31, 32, 33) along with Friedmann equations is shown below 



Redshift - Z 



FIG. 1; Higgs held as a function of redshift. 



Redshift - Z 



FIG. 2: Equation of state oJmggs for Higgs helds, as seen it starts with —1 then evolves 
towards quintessence regime after large enough time it comes back at —1. 

During the initial stages Z 3750, the evolution of the Higgs helds 772 , X 2 and are 
frozen by the expansion shown in hg. ([^ and acts as a negligibly small vacuum energy 
component with u = —1. As time proceeds, Higgs helds begins to evolve towards minimum 
of potential, the energy density in the Higgs helds start to dominate cosmologically (on 
Hubble scale). During the evolution equation of state parameter ujHiggs starts to increase 
and becomes > — 1 as shown in hg. ([^. In the very late (future) Universe {Z < 0), the helds 
comes to rest at the minimum of the potential and a period with a; = — 1 is reachieved to 
give an exponentially accelerating Universe. Since, oJmggs ~l/3 at any time in evolution, 
after cueff gets less than —1/3 in the evolution, an ever accelerating Universe is obtained in 
this model. 

As discussed before, Higgs held stability is needed to get an ever accelerating expanding 
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Universe. This was obtained by imposing the discrete Z 2 symmetry. The lightest Higgs helds 
r ]2 and X 2 in this case do not decay into any other Higgs held (since these helds are lighter 
than SM like and charged Higgs) and also into fermions (since they do not couple to them). 
Thus, we have a model of ever accelerating expanding Universe. 


0.4 
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Sc -0.2 
^"-0.4 
- 0.6 
- 0.8 
- 1.0 


FIG. 


The cjeff in the hg. (|^ starts from 0.167 (set by initial conditions ^Higgsi^t = 0 and 
^NRi„t = ^Rint = 0-h; NR: non-relativistic matter and R: relativistic matter) and decreases 
as the relativistic matter’s energy density decreases. This has been shown in hg. g for UJeS 
and in hg. (|^ for energy densities and comes down at 0. This and before is the time period 
when non-relativistic matter dominates and expanding Universe decelerates with highest 
rate {^Inr ~ 1) as non-relativistic matter domination pulls every thing inwards more than 
outwards Higgs negative pressure. After that Wes starts to decrease as the non-relativistic 
energy density decreases and Higgs relic energy density increases as shown in hg. 0 and 
([^ this time and afterwards Higgs negative pressure starts to dominates for forever (because 
second phase transition never occurs in our model) and Wes eventually settles down at —1. 
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3: Ehective equation of state Wg// = Qde^de + ^r^r + ^m^m- 
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FIG. 4: Energy densities of different components as a function of redshift. 
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FIG. 5: Energy densities of different components as a function of redshift. 


Figure (|^ and ([^ show the evolution of different components of the Universe as a function 
of redshift Z. In Fig. Q we can see the evolution of different components from Z 3750 
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to Z — 1 while fig. ([^ shows the evolution from Z = 1 to Z ^ —1. From fig. ([^, we 
note that the energy density decreases as a function of a~^ for non-relativistic matter 

(relativistic matter) while from fig. (|^ it can be seen that for high redshifts the energy 
density for Higgs fields remains approximately constant. However, for low redshifts (from 
fig- ([^) the Higgs energy density decreases to settle down at their minimum value {Eyac 
given by eq. (|^) as the fields come to rest at low redshifts (shown in fig. 0 ). 
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FIG. 6: Relic densities of different components as a function of /n[ao/a]. 
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FIG. 7: Relic densities of different components as a function of redshift. 


Figure (|^ is plotted for the evolution of the relic density of the non-relativistic matter, 
relativistic matter and the dark energy Higgs fields against ln{ao/a{t)), starting from the 
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non-relativistic and relativistic matter equality; where t is time. The hgure shows that only 
now the dark energy Higgs has started to dominate in the critical density while in the past 
non-relativistic matter was dominating. Figure ([^ is plotted for the evolution of the relic 
density of the non-relativistic matter, relativistic matter and the dark energy Higgs helds 
against the redshift Z starting from the non-relativistic and relativistic matter equality to 
Z ^ —1] where t is time. The hgure shows that for Z < 6 the dark energy Higgs has started 
to become prominent. At Z ^ 0.3235 the non-relativistic matter and relativistic matter 
energy densities are equal. For Z > 0.3235 non-relativistic matter is dominating and for 
Z < 0.3235 dark energy Higgs started to dominate causing the accelerated expansion of the 
Universe. 


It is quite clear from the plots that component Higgs helds of second Higgs doublet are a 
possible candidate for the dark energy if the present Universe happens to be evolving in the 


vacuum states described by inert doublet model given by the vacuum solution eq. (22). 

One needs to know that the initial conditions for charged held were taken in accordance with 
the observation that dark energy (vacuum) is not charged. 


V. DECAY(S) OF HIGGS FIELD(S) 


At tree level the doublet 02 in inert doublet model discussed above whose component 
helds are the candidate for dark energy is only coupled to the doublet <pi. The doublet 0i 
acts in an identical way as SM Higgs doublet. 

The interaction Lagrangian of Higgs helds of doublet 02 with Higgs helds of doublet 0i and 
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gauge bosons (extracted from eq. (§) is 


+ A 4 - \5)X2^V1 + ^^i(A3 + A 4 + \5)V2^V1 + A3Z^1^102 02 
+ ^(A3 + A 4 — A5)x2^^1^ + ^(A3 + A 4 + A5)772^^1^ + 2^3Vl‘^4>2 4’2 




(92^+91)"" "“■■'• 
^2^2 (^2 - ^2^) 
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2 , 9',V-2 S^2 


4(<7^2 + gl) 


( .2 ^ ,) 4^2"91„Z„ + jW-,W+,{r,i + xi + 


2„r 


2 22 / 


+ 0 n 'hM'l’tw-, + .>2-11-'%) + i-^S^=x29i,(.#.2-ir+„ - <t>tw-,) 

^^{ 92 + 92 ) 2 ^(^ 2 +^ 5 ) 

/ 2 / 2 
- ^ + 02“^+ m ) - ^ 2. Y 2 ZM 2 W+, - 02+^“ m )- 


2 \/ (5'2^ + 92 ) 


2 \/ (5'2^ + 92 ) 


( 35 ) 


To suppress the interaction of Higgs fields 772, X2 and 0 ^ with the gauge bosons the idea is 
that the SU( 2 ) doublet 02 is very weakly (and different than 0 i) coupled with the gauge 
bosons, thus here g2 7?! and g2 g'l- Thus, the decay modes that include the Higgs 
helds of the second Higgs SU( 2 ) doublet are negligible as compared to the modes that include 
the SM Higgs. 

The Higgs decay to a pair of (Higgs) scalar while considering only the on-shell width, the 
decay width is given by [ 19 ] 


T{H, HjHk) = (2 - 6,k)mH, 
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CuiHjHk is the coupling of different Higgs bosons Hi, Hj and 


From eq. ( 35 ) C, 
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The decay rate of gi to 772772, X2X2, 02 ^2 masses used in 
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the cosmological evolution determination in the previous section are, 


nm - 

-t X 2 X 2 ) 

= 1.53794 X 10-6 GeV 


V 2 V 2 ) 

= 1.53794 X 10-6 GeV 

r(>)i - 

0^ 02 ) 

= 0. GeV 


Also the decay rates for SM like Higgs boson is given in the graph below as a function of 
masses m,, 2 , and 
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From the above plot, we see that the decay rate for SM like Higgs to other Higgs pairs is 
peaked at mx = 51.3168GeV. For mx < 51.3168GeV the decay rate increases from zero to 
F = 15.4557GeV and F = 7.72786GeV for charged pair and neutral pair respectively, after 
that it decreases sharply and goes to zero at mx = 62.85GeV. 

It should also be mentioned that total decay width of SM like Higgs boson here is well 
within bounds of mass resolution 12 x 10 “^GeV of LHG [2S] (this is only the current results 
and the uncertainty in the measurements is high [2S] which is expected to be rehned [could 
come to SM expected value as more data is collected]) for the masses used in the cosmological 
evolution of the Higgs helds. One should also mention that the SM prediction of total decay 
width for Higgs boson is 4.21 x 10“^GeV with mass 126GeV [27] and is 4.07 x 10“^GeV with 
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FIG. 8 : Decay rates for rji —> xx 
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mass 125GeV |26] . In this model, we get three more decay channels of SM like Higgs which 
are to the other Higgs bosons pair. 

It should also be mentioned that the dark energy density does not modify as a consequence 
of the SM-like Higgs decay into dark energy Higgs pair for the masses of the Higgs fields taken 
here. Had we been able to achieve accelerated expansion with the mass of the second Higgs 
> 125GeV (which we could not because of fine tuning problem) the dark energy density 
would change significantly at redshift corresponding to the equilibrium temperature (freeze 
out temperature of dark energy Higgs) Teg = — 2mHsM- ^^at case, we could have 

a Universe whose different regions would accelerate, decelerate or contract differently — a 
case usually called bubble Universes, multi-verses or parallel Universes. These decays do not 
imply any additional fine tuning. 


VI. CONCLUSIONS 


Scalar fields are among the possible candidate for the observed accelerated expansion 
of the Universe. In this article, I argued that the Particle Physics developed so far must 
have something in/beyond SM which will explain the observed accelerated expansion of the 
Universe and hence will serve as the dark energy candidate. This is also necessary for the 
unification of matter and vacuum energy. This unification is necessary in my opinion. This 
article took specifically this approach that dark energy actually is some form of scalar field(s) 
which is(are) present in the so called inert doublet model. 

I have found that if the present Universe is described by the vacuum given by an inert 
doublet model then the component scalar fields of the inert doublet 02 can be one possible 
candidate for the dark energy. Since the present contribution of the dark energy to the 
critical energy density is about 0.7, this value is obtained by taking the masses of GP-even 
field(s) very low. The most important thing is that with the initial conditions set the mass 
of the charged (0^) field become arbitrary if we release the parametric constraint given by 


eq. (29) which is obtained using phase transitions bound, this model will fit for any value 


of mass of 02. One needs to keep in mind that the values of masses were chosen arbitrary 
just to get dark energy relic density ~ 0.7, changing the values of masses mildly does not 
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changes the relic density mnch. We conld also have some other masses set which wonld give 
the same evolntion as shown below for above masses set bnt there is some bonnd on the 
parameters of the potential to get high damping in the helds (which restricts the masses to 
remain in certain order of energy). 

As discussed above even though the masses were chosen arbitrary to get flmggs ~ 0.7, to 
avoid the oscillations in the helds we have to take A3 + A4 + A5 ~ 0(10“^^^) and A3 + A4 —A5 ~ 
0(10“^^®) and thus the masses of second Higgs doublet very small to ensure that we get 
cueff = —1/3 only once in the history of the Universe. The obtained results suggests to probe 
the low energy behavior at the Universe’s current temperature (T = 2.73°K) of the standard 
model and its extensions to obtain new information and possibility of dark energy candidacy 
in Particle Physics. 

To suppress the interaction of the Higgs helds 772, X2 and (pf with the gauge bosons Z 
and A{fj,), we imposed the condition that the SU(2) doublet 02 is very weakly (and diherent 
than 01) coupled with the gauge bosons. Thus here g 2 gi and 5^2 ^ 9 i- This allowed us 
to conclude that the decay modes that include the Higgs helds 772 , X 2 and pf are negligible 
as compared to the modes that include the SM Higgs. The couplings gi and g[ have the 
values of the SM Higgs gauge couplings in the SM but the values of g 2 and 5^2 niust be set 
using the results from LHC or any other accelerator. 

One thing that remains important to check in all extensions of the SM is whether the 
Higgs potential contains the vacuum instability or not? If it does, then does it make the 
vacuum instability more worse as compared to the SM or not? The answer to the question 
for our model is that although it contains the vacuum instability, due to the coupling of the 
second Higgs with the SM Higgs is (P(10“^^®), it will not affect the RGBs running of the SM 
Higgs. Thus will not make the vacuum instability worse. We expect the vacuum instability 
to occur at approximately the same scale as it occurs in the SM. 

When T becomes less than Tc^, symmetry in both helds is broken and all mixed helds of 
both doublets get massive (except Goldstone bosons). In this situation, one neutral Higgs 
can be made to acts in the same way as SM Higgs does (i.e. giving masses to other particles 
and having mass 125.7GeV) while other helds can be made to act as dark energy held(s) by 
setting the Yukawa interaction to cancel the other helds ehect on fermions. This should be 
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investigated in future. 

When describing a model for accelerated expansion of the Universe, it becomes essential 
to compare it with the ACDM model. In comparing our model with the ACDM, we note 
that on the cosmological scale our model is completely different from the standard ACDM 
and also very predictive. In our model ojHiggs is not constant over redshift which has been 
shown in £g. ([^. From fig. Q, we see that the model gives oJuiggs ~ —0.858 at Z = 0. 
Thus our model can be distinguished from the ACDM via the variation of OJHiggs from — 1 
on the cosmological scale. The deviation of oJuiggs at low redshift 2 < Z < 0 is given in 
fig- ([^. In the plot, we see that while we started from oJniggs = —1 at Z 3750, around 
Z 6 00 Higgs start to get greater than —1. At redshift Z = 0 the deviation of ooHiggs from — 1 
we obtain in our model for explicit values of parameters using MATHEMATICA software 
is 0.1419828 ± 3.2433649325 x 10“®. In concluding, since we get ooe/j < —1/3 after solving 
the Euler Lagrange equations numerically, the proposed Higgs held could cause the current 
observed accelerated expansion. 

It is also worth concluding that because of the type of Yukawa interaction the second 
phase transition (the case of m |2 ^ 0) will not change fermions phenomenology but due to 
the SU(2)(8)U(1) gauge invariance gauge bosons phenomenology will get changed (in gen¬ 
eral). Setting the coupling constants g 2 and g 2 might end up not changing gauge bosons 
phenomenology subsequently. This should be investigated in future theoretically. 
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